Momentum distribution functions in ensembles: the inequivalence of microcannonical 
and canonical ensembles in a finite ultracold system 
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It is demonstrated that in many thermodynamic textbooks the equivalence of the different en- 
sembles is achieved in the thermodynamic limit. In this present work we remark the inequivalence 
of microcannonical and canonical ensembles in a finite ultracold system at low energies. We calcu- 
late the microcanonical momentum distribution function (MDF) in a system of identical fermions 
(bosons). We find that, the microcanonical MDF deviates from the canonical one, which is the 
Fermi-Dirac (Bose-Einstein) function, in a finite system at low energies where the single-particle 
density of states and its inverse are finite. 
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I. INTRODUCTION 

The recent progress in ultracold atomic gases experi- 
ments provides an important test bed for studying the 
isolated quantum many-body systems, where these sys- 
tems are sufficiently weakly coupled to the external envi- 
ronment*. The non-equilibrium dynamics after quench- 
ing in quantum systems is one of the hot topics. In an 
isolated system of conserved finite particles, one may ask 
in what conditions the microcanonical and cannonical en- 
sembles are equivalent? 

According to statistical mechanics, the momentum dis- 
tribution function (MDF) in a canonical ensemble is 
the Fermi-Dirac (Bose-Einstein) function for a system 
of identical fermions (bosons )2. The Fermi-Dirac (Bose- 
Einstein) function is first derived in an ideal gas of 
fermions (bosons) , but now they are generally used in in- 
teracting systems in terms of Fermi liquid theory^, where 
the free fermions (bosons) are replaced by the emergent 
quasiparticlcs with rcnormalizcd energies. However, the 
MDF in a microcanonical ensemble is difficult to calcu- 
late directly. A textbook derivation of the Fermi-Dirac 
(Bose-Einstein) function in microcanonical ensembles^ is 
based on an assumption that the MDF is the most proba- 
ble distribution guaranteed by the principle of the equiva- 
lence of ensembles, which states that the macroscopic ob- 
servables (e.g., the MDF) in canonical and microcanoni- 
cal ensembles are the same in the absence of long-range 
interactions in the thermodynamic limit£~— . However, it 
is obliged to study in which conditions the equivalence of 
ensembles breaks down in a finite system. 

Recently, the MDF in cold-atom systems after a quan- 
tum quench causes great interests both in theory and 
experiments^—. In a system far from intcgrability, ac- 
cording to the eigenstate thcrmalization hypcrthesisi&~— , 
the steady state after a quench can be described by an 
ensemble of eigenstate a 18 ' , close to each other in energy. 
In other words, a non-integrable system will thermalize 
into a microcanonical ensemble^. At the same time, a 
typical cold-atom system in experiments contains only a 



few thousand particle s 10 ' 20 , far from the thermodynamic 
limit. The MDF of microcanonical ensembles in a finite 
system is now accessible in experiments implemented on 
ultracold atoms. Then it is necessary to revisit the MDF 
in microcanonical ensembles when the equivalence of en- 
sembles breaks down in a finite system. 

In a finite one-dimensional system, Rigol studied the 
inequivalence of the canonical and grand-canonical en- 
sembles 21 . In a system of non-interacting fermions a nu- 
merical comparison between the MDFs in microcanonical 
and canonical ensembles with equidistant single-particle 
levels is discussed at some special energies^. However, 
a systematic study of the inequivalence between micro- 
canonical and canonical ensembles is still lack. A par- 
ticularly interesting case is the low-energy state, when 
the system energy is comparable with the inverse of the 
single-particle density of state which is nonzero in a fi- 
nite system. Obviously, the MDF must be zero when the 
single-particle energy exceeds the system energy, indicat- 
ing that it is not a smooth function, different from the 
Fermi-Dirac (Boso-Einstein) function at finite tempera- 
ture. 

In this paper, we discuss the MDF in microcanonical 
ensembles. We suppose that the system is finite by set- 
ting a finite single-particle density of state. We calculate 
the microcanonical MDF (MMDF) by relating it to the 
canonical MDF (CMDF). A well-known relation between 
canonical and microcanonical ensembles is that canonical 
ensembles describe the subsystems of a microcanonical 
ensemble^. In this paper, we show that the MMDF of 
a system can also be calculated according to its CMDF. 
We find that the MMDF deviates from the CMDF at 
low energies, while becoming equivalent at high energies. 
Our results show the breakdown of the equivalence prin- 
ciple, at the same time provide a way to understand the 
MDF in cold-atom systems after a quench, which is be- 
lieved to be the MMDF, while difficult to address directly 
due to the lack of reliable analytical or numerical meth- 
dslU2i2tfZ. 

The contents of the paper are arranged as follows. In 
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Sect. [IT] we introduce the methods to derive the MMDF 
in a finite system. In Sect. IIIII the MMDF for fermions 
are calculated and the first- and second-order approxi- 
mations for the MDFs are derived. In Sect. II VI we report 
the corresponding results for a system of bosons. At last, 
a concluding section summarizes our results. 



II. METHOD 

Let us suppose a system of identical fermions (bosons) 
in a canonical ensemble. Its MDF is written as 
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(1) 



where /? = l/fc^T is the inverse of the temperature and 
e k the quasiparticlc energy. The plus and minus signs 
(±) are for fermions and bosons, respectively. 

The MDF is the expectation value of the momentum 
operator h k , and can be expressed as 
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where n k {ipi) = {4>i\n k \ipi) and \ipi) is the eigenstate 
with energy Ei. The partition function is defined as 
Z = ^i&^^ Ei - The sum is over all the eigenstates. 
When dividing the range of eigenenergies into small in- 
tervals [E — ^-,E + the sum over eigenstates in 
Eq. @ is rewritten as 



n k = 



E E 



(3) 



E Ei e[E-££-,E+±2-] 



where the sum to E is over all the energy intervals. Since 
AE is taken to be small, the eigenenergies Ei in an in- 
terval can be treated as a constant. Thus, the MDF is 
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where fl(E) denotes the number of eigenstates in the in- 
terval [E 
comes Z - 



AE 
2 



E + and the partition function be- 
J2 E e~ PE - The term in the bracket of Eq. (gj) 
is exactly the average of the microcanonical ensemble of 
the MDF, denoted as n k (E) mc . 

The definition of microcanonical ensemble demands 
that the number of eigenstates in the energy interval must 
be large enough to make the function n k (E) mc smooth. 
This condition can be naturally satisfied by converting 
the sum to E into an integral over E ranging from to 
oo. We set the ground energy to zero and suppose that 
low energy excitations in the system are gapless. The 
Eq. ([4]) becomes 



n k 



- jf dED(E)e-? E n k (E) mc , 



(5) 



where D{E) is the density of eigenstates and the partition 
function becomes Z(f3) = J Q dED(E)e~ l3E . Converting 
the summation into an integral smears out the fluctua- 
tion of MDFs in the eigenstates. This is equivalent to 
calculating the average with respect to a microcanonical 
ensemble. In a system far from integrability, the eigen- 
state thcrmalization hyperthesis can be applied^, then 
the f luctuat ion of MDFs is absent. In this case, the func- 
tion n k (E) mc is equal to the MDF of an eigenstate with 
energy E. The MMDF and the MDF of the eigenstate 
in systems that the eigenstate thcrmalization hyperthe- 
sis can be applied will not be distinguished next, both 
denoted by n k (E) = n k (E) mc . 

Comparing Eq. ([5]) with Eq. (JlJ, we finally obtain 



/•OO 

/ dED(E)e' pE n k (E) = 
Jo 



Z(f3) 

,0e k ± I ■ 



(0) 



This equation shows that the Laplace transformation of 
D(E)n k (E) is the product of the Fcrmi-Dirac (Bosc- 
Einstcin) function and the partition function determined 
by D(E). The MDF n k (E) can then be obtained by an 
inverse Laplace transformation. The MMDF is a univer- 
sal function of energy, depending only upon the density 
of eigenstates in the system. 

The equivalence of ensembles can be understood as 
a result of that the product D{E)e~^ E should be a 5- 
like function in the thermodynamic limit. In this paper, 
however, we consider a finite system, where the single- 
particle density of states, proportional to the volume of 
the system, is finite. Then we choose the inverse of the 
single-particle density of states as the energy unit. We 
find that the equivalence of ensembles breaks down at low 
energies by showing that the MMDF is different from the 
Fcrmi-Dirac (Bose-Einstein) function. 

Next we discuss the MMDF in fermionic and bosonic 
systems, respectively 



III. MMDF IN A SYSTEM OF FERMIONS 

We first discuss a system of identical fermions. The 
density of eigenstates D(E) is calculated. In the ground 
state, all the quasiparticle levels lower than the Fermi en- 
ergy (set to zero) are occupied and the others are empty. 
The excited eigenstates are classified according to the 
number of the quasiparticlcs excited. The total density 
of eigenstates is expressed as an infinite series, 



D(E) = J2D n (E), 



(7) 



n=l 



where D n (E) denotes the density of eigenstates in which 
?i particles are excited. Then D n (E) is expressed as 



D n (E) 



dE 
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FIG. 1. The MMDF in the lowest order approximation in a 
system of fermions. Here E denotes the difference between 
the system energy and the ground state energy, and Eh the 
quasiparticle energy. The Fermi energy, defined as the energy 
of the highest occupied quasiparticle state in the ground state, 
is set to zero. 



where 0(efe) is the quasiparticle density of states and 
8(e) denotes the Heaviside function. We concentrate on 
the low-energy eigenstates, in which only the quasiparti- 
cles close to the Fermi level contribute to the function 
D n (E). In the wide-band limi^, ()(e k ) close to the 
Fermi level is a constant denoted by T>. Then we find 

Only the 



(2n-l)! ■ 

The function D(E) is a power series of E. 



first several terms are important in calculating the MDF 
at low energies. We first take the lowest order approxima- 
tion D(E) — T> 2 E, and find that the partition function is 
Z(fi~) — correspondingly. The lowest order approxi- 
mation is valid in the limit /3 — ► oo or E — > 0. By solving 
Eq. we find that the corresponding MMDF is a step 
function: 



n(e k ,E) 



if 



2j+2 



2j+l 



with j = 0,1,2, and rescaled dimcnsionless energy 
e k = e k /E > . The MMDF as e fe < is given by 
n(e k ,E) = 1 - n(-e k ,E). In Fig. 1, we plot n(e k ,E) as 
a function of e k . The MMDF is very different from the 
CMDF which is a Fermi-Dirac function. It shows quan- 
tized features. We find n(e kl E) = 1 for e k < —E, and 
n(e k ,E) = for e k > E, respectively. This truncation 
is natural because the quasiparticle excitation with an 
energy larger than the eigenstate itself is forbidden. In 
the regime — E < e k < E, the MMDF has a sawtooth 
shape and becomes irregular close to the Fermi energy 
due to the quantization of quasiparticle excitations in a 
finite system in which the quasiparticle energies are dis- 
crete. In the lowest order, the MMDF depends only upon 
the ratio e k /E. It recovers the well-known ground state 
MDF 6 (-e fc ) in the limit E -> 0. 

We then consider the second order term in D(E) and 
have D(E) = V 2 E + Now the partition function 

is Z = ^ + ^. The MMDF as e k > is found to be 



ri ■,. |1 , 1+P 2 [1-(J- + l)hf /6 + j(j + l)V 2 ~ej/2 

[1 " (k{] + 1)1 l+W/6 ' 

- . 1 + V 2 [1 - (j + l)gfe] 2 /2 + V 2 h [1 - (j + l)g fc ] /2 + f&ej/e 

£fe 7 = 

l + P 2 /6 



n{e k ,E) 



if ^ < h < 



2j+2 



2 j + 1 



' if 27TT ^ gfe < h 



, (10) 



where T> = T>E is a dimensionless quantity. The n(t kl E) 
as E > depends not only on e k but also on V after 
considering the second order term in the density of eigen- 
states. In Eq. (|10p . the system energy E is contained in 
T>, which increases linearly with E. And Eq. (|10p returns 
to Eq. © as T> — > 0. Since D(E) is a power series of T>, 
the second order approximation is valid as T> is not too 
large, i.e., E is not much larger than 1/2?. This defines 
the regime of system energies in which Eq. (|10[) is cor- 
rect. In Fig. 2, we plot the MMDF in Eq. KTDJ} at different 
E. As the energy increasing, the sawtooth shape in the 
MDF is softened. The MMDF gradually changes into the 
Fermi-Dirac function, i.e., the CMDF. In fact, the differ- 
ence between the MMDF and CMDF at high energies will 
totally disappear after considering higher order terms in 
D(E), a result of the equivalence of ensembles. We take 



the large E limit (V — > oo) in Eq. (fT0|) . and then take 
the limit e k -> 0, and find n(e k ,E) = | — ji k + 0(e|). 
The slope of the MMDF at the Fermi energy is — 3/(4E). 
We define an effective temperature so that the slope of 
the corresponding Fermi-Dirac function at the Fermi en- 
ergy is — 3/(4E). The effective temperature is found to 
be T = £7(3fcs). Because Eq. (fl"U| is valid at low en- 
ergies, this effective temperature should be attached to 
microcanonical ensembles of a system of fermions at the 
low-energy limit. 

The MDF in microcanonical and canonical ensembles 
are very different at low temperatures. We have ex- 
pressed D(E) in a power series of E. Correspondingly, 
n(e k , E) is a power series of T>. The result in Eq. (fTU|) is 
kept to second order of T>. As T> is small, the 0(T> 2 ) 
terms in Eq. JTDJl is negligible. The MMDF has the 
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FIG. 2. The MMDF at different energies in a system of 
fermions. Here E denotes the difference between the system 
energy and the ground state energy, and e k the quasiparticle 
energy. And Ek — denotes the Fermi energy. The canonical 
MDF at the effective temperature T = E/(3Icb) is plotted as 
a comparison, labeled by the acronym CMDF. 

sawtooth shape similar to that in Eq. ([9]), clearly distin- 
guished from the Fermi-Dirac function. At low energies, 
the quantization effect of quasiparticle excitations is im- 
portant, resulting the oscillation behavior of the MDF in 
the regime — E < e k < E. 

The steady MDF of a nonintegrable system driven out 
of equilibrium is believed to be the MMDF—. And the 
corresponding energy E is decided by the initial condi- 
tion. While, the CMDF is obtained in an open system 
contacting with a large thermal reservoir. In this pa- 
per, we find that the equivalence of ensembles is broken 
in a finite system at low energies E ~ 1/2?. The in- 
equivalence of ensembles has been intensively discussed 
in classical mechanics in the presence of long-range inter- 
actions^— iSSLSfi. Our results must be distinguished from 
those discussions, because the inequivalence of ensembles 
in our paper is a finite-size effect due to the quantization 
of excitations. The equivalence will be recovered at high 
temperatures. 

Our results are related to experiments on isolated ul- 
tracold atoms, in which the MMDF can be directly mea- 
sured. The MMDF in Eq. (fTU)) can be observed in a 
system in which the CMDF is the Fermi-Dirac function, 
e.g., in a Fermi gas or a Fermi liquid. This suggests 
that the system is two- or three-dimensional, since the 
one-dimensional fermions are described by a Luttingcr 



liquid^. Here we discuss the condition that the inequiv- 
alence of ensembles is obvious in two dimensions. The 
Fermi gas in two dimensions has been realized in ex- 
periments (see Refsi22r— for examples). We choose an 
isotropic dispersion relation e k — ^kBr, where m denotes 
the quasiparticle mass. In two dimensions, the quasipar- 
ticle density of states is a constant as V = j^p-, where 
S is the area of the system. Then we have T> — f-Sf. 
We use the relation between the temperature and the sys- 
tem energy obtained above and finally we get T = J^^g • 

The inequivalence of ensembles is obvious as d ~ 1 . Then 
the temperature must be inversely proportional to the 
system area. 

When considering a system of 40 K atoms^ with the 
mass 6.6 x 10~ 26 kg and a typical area of two-dimensional 
Fermi gas S = l^im 2 , we obtain the typical temperature 
T ~ 2.5 x 10~ 8 K where the inequivalence of ensembles 
can be observed. 



IV. MMDF IN A SYSTEM OF BOSONS 

Next we discuss a system of bosons. Without the Pauli 
exclusion principle, the infinite number of bosons occu- 
pies the lowest quasiparticle level (set to zero) in the 
ground state. Then the density of excited eigenstates 
becomes 




Here we consider the case that the quasiparticle density 
of states is a constant 2?, e.g., in two-dimensional sys- 
tems. Then the density of eigenstates is found to be 

D ( E ) = J2n=i V (n-"iy. ■ A g ail h tne function n(e k ,E) is 
a power series of T> = VE and can be calculated order 
by order at low energies. We keep D(E) to order 0(E 2 ): 
D(E) = V + V 2 E + V 3 E 2 /2. This corresponds to a func- 
tion n(e k , E) kept to order V 2 . The partition function is 
Z = ? + + jp. The Laplace transformation in the 
case of bosons is defined as 

J™ dEe-P E D{E)n{e k ,E) = -jg^L (12) 
Performing the inverse Laplace transformation, we find 



1+2? 



n{e k ,E) = j- 



i±Lf)~e k + \ ( (1 - j~e k ) 2 + (J - l)h(l - jh) + 



(j-l)(2j-l) - 2 



1 + V + V 2 /2 

I 



(13) 



as < i k = e k /E < \ with j = 0, 1, 2, ■ 



In Fig. 3, we plot the bosonic MMDF at different ener- 
gies. The MMDF is a step function at low energies, while 
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FIG. 3. The MMDF in a system of bosons at different en- 
ergies. Here E denotes the difference between the system 
energy and the ground state energy, and ek the quasiparticle 
energy. And eu = denotes the energy of the lowest quasi- 
particle state. The function n(ck,E) = — | + ^ in the 

limit T> — > oo is plotted as a comparison, labeled by 'limit'. 



power law decay is that T> is not too large so that the 
second order approximation to n(ek,E) is effective. 



V. CONCLUSIONS 

In summary, we calculate MMDFs in a system of 
fermions and bosons. The MMDF given in this paper 
is universal in systems where the corresponding CMDF 
of quasiparticles is the Fermi-Dirac (Bosc-Einstcin) func- 
tion and the low energy excitations are gapless. The 
MMDF is found to be different from the CMDF at low 
energies as E ~ l/T>, where V is the quasiparticle density 
of states. It shows a sawtooth shape for fermions and a 
power-law decay for bosons. Our results show the break- 
down of the equivalence of ensembles, as a result of the 
quantization of quasiparticle excitations which can only 
be observed in small systems, since T> is proportional to 
the system volume. The MMDF obtained in this pa- 
per is expected to be verified in experiments of isolated 
cold- atom systems, which is believed to thermalize into 
a microcanonical state after a quench. 



changes into a smooth one as the energy increasing. Tak- 
ing the limit V — > oo, we get n(ek,E) = gi- — \ + ^. 
This indicates that at low energies the bosonic MMDF 
as a function of quasiparticle energy decreases in a power 
way instead of in an exponential way suggested by the 
Bose-Einstein function. The condition of observing this 
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